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Abstract 



Work of the last few years has shown that the key algebraic features of Logarithmic Con- 
formal Field Theories (LCFTs) are already present in some finite lattice systems (such as the 
a ' XXZ spin- 1/2 chain) before the continuum limit is taken. This has provided a very convenient 

way to analyze the structure of indecomposable Virasoro modules and to obtain fusion rules 
for a variety of models such as (boundary) percolation etc. 

LCFTs allow for additional quantum numbers describing the fine structure of the indecom- 
■ posable modules, and generalizing the 'fe-number' introduced initially by Gurarie for the c — 

case. The determination of these indecomposability parameters (or logarithmic couplings) has 
given rise to a lot of algebraic work, but their physical meaning has remained somewhat elusive. 
In a recent paper, a way to measure b for boundary percolation and polymers was proposed. 
I We generalize this work here by devising a general strategy to compute matrix elements of 

Virasoro generators from the numerical analysis of lattice models and their continuum limit. 
The method is applied to XXZ spin- 1/2 and spin-1 chains with open (free) boundary condi- 
tions. They are related to gl(n -I- m\m) and osp(n + 2m|2m)-invariant superspin chains and to 
non-linear sigma models with supercoset target spaces. These models can also be formulated 
in terms of dense and dilute loop gas. 

We check the method in many cases where the results were already known analytically. 
Kji • Furthermore, we also confront our findings with a construction generalizing Gurarie's, where 

' logarithms emerge naturally in operator product expansions to compensate for apparently 

, divergent terms. This argument actually allows us to compute indecomposability parameters 

in any logarithmic theory. A central result of our study is the construction of a Kac table for the 
indecomposability parameters of the logarithmic minimal models £A4{l,p) and £,A4{p,p+ 1). 



Introduction 

Logarithmic Conformal Field Theories (LCFTs) are characterized by the fact that the dilatation 
operator Lq is non-diagonalizable and has a Jordan cell structure. This property leads to the ap- 
pearance of logarithms in correlation functions and to indecomposable operator product expansions. 
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Such theories have emerged recently in many physical situations such as geometrical problems, dis- 
ordered free fermions models in 2+1 dimensions or the AdS/CFT correspondence in string theory. 

Geometrical applications include self-avoiding random walks (polymers) or percolation (see e.g. 
[1112]). Although it would be fair to say to these problems are well understood, a fully consistent 
field theory describing geometrical observables {e.g., four point correlation functions in the bulk) is 
still missing. Another very exciting question, with considerable applications to condensed matter, 
concerns the critical point in non-interacting disordered fermion models in 2+1 dimensions, which 
are believed to be described by c = LCFTs. Physical realizations of such systems are given 
by transition between plateaux in the Integer/Spin Quantum Hall Effects. The appearance of 
logarithmic correlators in such theories cannot be avoided [3] and a deep understanding of LCFTs 
is necessary in order to describe their low-energy physics. While much has been done about the 
Spin Quantum Hall Effect, as it can be related to the classical percolation problem [1], very little is 
known about the field theory describing the transition in the Integer Quantum Hall Effect (see |S] 
for a review). Logarithms in CFTs also arise naturally in the context of non- linear sigma models 
with super target space. These quantum field theories play a major role in the AdS/CFT duality. 
For example, the PSL(2|2) sigma model is related |5] to strings living in AdSa x S^. 

The key feature of logarithmic CFTs, as opposed to simpler non- unitary theories (such as the 
Yang-Lee singularity, or other non- unitary minimal models), is indecomposability. This property 
was probably observed first by Rozansky and Saleur [7] who studied a WZW model with U(l|l) 
supergroup symmetry and vanishing central charge. They related the non-simplicity of U(l|l) to 
the possibility of non-diagonalizability of Lq and logarithmic dependence in four-point correlation 
functions. The study of logarithmic c = theories in a more systematic fashion then begun with 
a serie of papers by Gurarie and Ludwig. Gurarie [S] first noticed that logarithmic operators were 
necessary in order to construct a consistent field theory at c = 0. A similar observation was made 
by Cardy using a replica approach [3]. Gurarie [2] and Gurarie and Ludwig [TU] then related 
the existence of a logarithmic partner t{z) for the stress energy tensor T{z) to ill-defined terms 
in operator product expansions. They introduced one of the firs10 indecomposability parameters, 
usually denoted &, which at the time was interpreted as a new 'anomaly' that would play the role 
of a central charge when c = 0. Using some heuristic arguments, they were also able to predict two 
possible values, ^ = — § smd 6 = |, which would distinguish between two fundamentally different 
LCFTs. 

Instead of a single parameter b, it is now well accepted that a LCFT is characterized rather by a 
complex structure of indecomposable Virasoro modules, with a infinite number of indecomposability 
parameters needed to describe the whole pattern. Two lines of thought have been considered. The 
first one is to deal directly with abstract indecomposable Virasoro modules to try to understand 
and classify their structures [HI [T^l [21 1131 E] • Progress has been steady in this direction, and the 
module involving T and its partner t at c = now appears as a particular case. These algebraic 
studies have led to new predictions for indecomposability parameters for various logarithmic pairs. 
The second idea is somehow more concrete, and consists in studying directly lattice models which 
can be thought of as lattice regularizations of LCFTs. For specific values of the parameters, the 
lattice Hamiltonians are non-diagonalizable and have a Jordan cell structure that mimics that of the 
continuum theory. This was mainly done by Read and Saleur and independently by Pearce, 
Rasmussen and Zuber [THl HZl [T5] . This approach essentially relies on the 'similarity' between the 
indecomposable modules of the Temperley-Lieb algebra and those of the Virasoro algebra. Both 

^Several other indecomposability parameters were computed by Kausch and Gaberdiel a few years earlier 1121 . in 
a different context. 
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approaches yield a consistent picture of boundary (chiral) LCFTs, with common algebraic structure 
and fusion rules deduced from both methods. As for the similarity, it is now better understood in 
terms of common quantum group structures [TO] , 

While the global structure of the Virasoro modules in the case of boundary LCFTs is slowly get- 
ting under control, many questions remain about the existence and values of the parameters arising 
from indecomposability. It is, moreover, tempting to think of these parameters as extensions of the 
structure constants of the operator product algebra of the theory, and thus to wonder what they 
encode physically, and whether it is possible to access them, experimentally or at least numerically, 
in the context of lattice simulations. Recently, Dubail, Jacobsen and Saleur [SOj suggested a con- 
crete method to measure one indecomposability parameter: the ?>-number of Gurarie in the case of 
c = theories. Their method involved a c = specific approach, the so called 'trousers trick', and 
led to the measurement of b for percolation and polymers. The observed values — b = — | and b — ^ 
respectively, were found in agreement with the predictions of Mathieu and Ridout [5] . 

The method of [20j is geometrically very appealing, but unfortunately does not extend to many 
other LCFTs, nor does it allow the study of indecomposable modules occurring at larger values of 
the conformal weight. It also does not seem to be generalizable to the bulk (non-chiral) case. On 
the other hand, having at one's disposal a 'probe' to investigate the detailed structure of modules 
in LCFTs seems rather essential to make progress in this, so far, very abstract topic. We have 
thus reconsidered the problem using a different route, which involves the identification of a lattice 
stress energy tensor. The original idea for doing so goes back at least to a paper by Koo and Saleur 
PT| . who themselves generalized the pioneering work of Kadanoff and Ceva [22]. The upshot of 
the proposal in [5T] was that the Temperley-Lieb algebra does not only reproduce in finite size 
indecomposable modules that mimic exactly several indecomposables of the Virasoro algebra: on 
top of this, there exist (infinitely many) linear combinations of words in the Temperley-Lieb algebra 
whose action in these modules reproduces, once properly interpreted, the action of all the Virasoro 
generators in the thermodynamic hmit. It thus should be possible, in principle, to reconstruct from 
the lattice all states of the boundary LCFT, and to measure all the matrix elements of all the 
Virasoro generators, hence to determine 'experimentally' all the information about the Virasoro 
modules in the model. 

Of course, this program is very difficult to implement in practice. Thanks in part to the progress 
accomplished in [Ellin], it is however not impossible, and this is what we do in this paper. 

Once the general strategy is under control, it turns out that we can study many more cases than 
those considered so far in papers on indecomposable Virasoro modules [TH [T31 [H] . These authors 
used mainly two different methods to compute indecomposability parameters. One is based on the 
so-called Nahm-Gaberdiel-Kausch algorithm to compute fusion products between indecomposable 
Virasoro modules [12]; while the other, used in Refs. [EKH], consider different quotients of 'glueing' 
of two Verma modules, incomposability parameters are then obtained as a solution of singular- vector 
equations. Note also that they can be computed using SLE considerations [23], or even directly 
by solving the differential equations satisfied by several four-point correlation functions. The latter 
method was one of the tools used by Gurarie and Ludwig [21] to study the case of the stress 
energy tensor at c = 0. While it should be possible to generalize the methods of these standard 
references to analyze the Jordan cells and calculate the indecomposability parameters in our cases, 
we find it more convenient here to use a 'short cut' analysis based on a generalization of the original 
arguments of Gurarie for the c = case. This is discussed, as a preamble, in the first section. The 
general framework is then discussed in section 2, where we introduce the lattice models, and use and 
generalize the algebraic arguments of Ref. [E] to deduce general indecomposable structures. The 
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next two sections are then dedicated to numerical measurements of indecomposability parameters. 
All our results are consistent with the predictions in [T^ [T31 [H] , but we obtain many results beyond 
these references — those for the first few logarithmic minimal models £Ai{l,p) and £A4{p,p +1) 
are summarized in table El A few conclusions are finally gathered in the last section. 

1 Preamble: indecomposability parameters and e — > argu- 
ment 

In this section, we define in simple terms the indecomposability parameter associated with a general 
pair of logarithmic operators. We then extend the c — > argument of Gurarie to other logarithmic 
CFTs and we show that this allows us to predict the value of the indecomposability parameters of 
a given theory. 

1.1 Jordan cells and indecomposability parameters 

It is now well known that the appearance of logarithms in correlation function is related to the 
non-diagonalizability of the Lq operator ^8 . This Jordan cell structure of the Hamiltonian (Lq) is 
itself deeply related to the non-semisimplicity of the underlying symmetry algebra of the theory. 
We will come back to these algebraic considerations later and choose here to focus on the Jordan 
cell structure. Consider a pair of logarithmic operators {4i{z),ip(z)) with conformal weight h that 
are mixed by Lq into a Jordan cell. In the basis {(j>,ip), the generator of the scale transformation 
reads 

Invariance under global conformal transformations then fixes [H^ the form of the correlation functions 

(0(z)0(O)) = 

{Hz)m) = 
mz)^m = 

where 9 and /3 are two parameters. While the constant 9 is arbitrary and can be canceled by a 
choice ^ ij) — ■§p4>^ the parameter /3 is a fundamental number that characterizes the structure 
of the Jordan cell. It is also important to remark that must be a null-field by conformal 

invariance, that is to say, introducing the usual Virasoro bilinear form, {4)\(t>) = 0. Actually, we 
know more about the algebraic structure hidden behind the non-diagonalizability of Lq. As we will 
see in details in the following, the fields 'ip{z) and (/)(z) always appear at the top and the bottom of 
a larger structure called a staggered Virasoro module [51[TTJ[T3], also called projectiv^ module in 

^Note that this nomenclature is somewhat dangerous as strictly speaking, these modules may not be projective 
over the Virasoro algebra in the mathematical sense. However, they can be seen as 'scaling limit' of projective 
modules over the Temperley-Lieb algebra that arises in lattice models (c/. next section), it is thus very tempting to 
call them projective anyway. 



y2h 



9-2f3\ogz 



(1.2a) 
(1.2b) 
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Ref [T3]. The fields in these modules are organized in a diamond-shaped structure that we note 



V= £, p . (1.3) 

\ ^ 

The arrows represent the action of Virasoro generators so that the whole structure can be induced 
by action of the L„'s on the field ^^{z) while 0(z) belongs to an invariant submodule. The whole 
module is reducible but indecomposable under the action of the Virasoro algebra. Remark that Lq 
maps -0 onto due to the Jordan cell action. The conformal weights of the different fields in p.3p 
satisfy < h ^ ^ h^i, < hp. The field (j>{z) is a null descendant of ^ 

<l>{z) = A^{z), A = L_„ + a(i)L_„+iL_i + --- + a(p(„)_i)L'Ii, (1.4) 

where n — h — and P{n) is the number of partitions of the integer n. The Q;(i) coefficients are 
uniquely fixed by the null-vector condition L^icj) — — 0. 

If two fields tp{z) and (/)(z) satisfy all these relations, we say that ip{z) is the logarithmic partner 
of the null-field (j){z), with indecompos ability parameter (also called logarithmic coupling) /?. Using 
the Virasoro scalar product, we see that 

p^{m- (1-5) 

where we normalized ^(z) such that (^|^) = 1. Note that it is also possible to define P through the 
equation 

AH{z) = Piiz). (1.6) 

It is important to notice at this point that the choice that we adopted for the normalization of the 
operator A is crucial for the value of /3. Different choices have been used in the literature, and 
some of them may yield simpler expressions for /3. Unless otherwise indicated, we always use the 
convention given by eq. I|1.4p . which is consistent with the normalization of the stress energy tensor 
in c = theories: T = L-2I so A = L_2 (see next paragraph). 

Finally, let us remark that the field p{z) never enters the computations of j3 so that it can be 
ignored as far as indecomposability parameters are concerned. 



1.2 c — 7- catastrophe and the stress energy tensor 

We now show that the /3's are fixed by a very simple argument relying on operator product expan- 
sions (OPEs) . The basic idea was developed in [SI [TD] , and a similar discussion may be found in 
. Consider a general CFT with the following parametrization of the central charge and of 
the Kac formula 

c = 1- , (1.7a) 

xyx + 1) 

_ [(x + l)r - xsf - 1 
^'^'^ " 4x(a; + l) • ^ 
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Conformal invariance fixes the OPE of an operator ^h{z) with itself to be of the form 



~/2h 



2h 



^r(o) 



(1- 



where T{z) = L-2I is the stress energy tensor of the theory. This expression is clearly ill-defined as 
c — > 0. In a nutshell, the idea of Gurarie was to introduce another field $1,5 for percolation (resp. 
$3^1 for dilute polymers) with conformal weight ht = /11.5 = 2 (resp. ht = ^13,1 = 2) at c = to 
cancel this divergence!!. Let us focus on the percolation case. When c is slightly different from 
= 2 + e), we can normalize^ the field $1.5 such that the OPE reads 



$,(z)$,(0) ~ ^ 



1 + — z2t(0) + z'''$i5(0) 
c 



We then define a new field t(z) as 



, , 2/1 (TIT) , , 2h^, , 
cp(e) c 



(1.9) 



(1.10) 



where /3(e) = --^^ and (T|T) 
defined as c — J' 



<i>,(z)$,(0) ^ ^ 



|. The OPE then involves quantities that are perfectly well 

^ r h ^ 

(1.11) 



l + ^z2(T(0)logz + <(0)) + 



with /3 = lime^.0 /3(e)- It is important to realize that the new field t{z) is perfectly well defined as 
c — >■ 0, while the Lq eigenvector $1^5 is not. In particular, one can then check (see e.g. [25|) that 
the fields T(z) and t{z) then satisfy the standard equations for logarithmic operators 



(T(z)T(O)) 
{T{zm) 

mm) 



^ 



0-2/31ogz 



(1.12a) 
(1.12b) 

(1.12c) 



with 6 a constant. A straightforward calculation using eq. (|1.10p also shows that Lot = 2i + T as 
expected. Note that it is quite general that one of the eigenvectors (here ^1,5) of an operator (here 
To) diverges as one tunes a parameter (here e) to approach an indecomposable point. One can then 
construct a new Jordan vector (here t) by canceling the diverging part in the ill-defined eigenvector 
by taking an appropriate combination with the eigenvector that has the same eigenvalue at the 
indecomposable point (here T). 

T(z) and its logarithmic partner t{z) are a special case of the general structure described in the 
previous paragraph. In particular, they are organized in a diamond structure like (|1.3p with ^ = /, 
(j)^T,^ = t, A = T_2 and L2t = /3T 



•^Actually, the original guess of Gurarie and Ludwig was #1,5 for dilute polymers and $3,1 for percolation. It shall 
become clear in the next sections that this result was correct up to a switch, as already remarked in Refs. [2l 1201 . 

^For the sake of clarity, we absorb the coefficient C^^^ /ag, into the normalization of ^i.s as this will play no 
role in the following "T^^r, ^^^ffi^i^n+o r'k 

(J* — <^<i, <!, • 



The coefficients C*^ are the usual structure constants that appear in OPEs; note also that 
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Gurarie and Ludwig then inferred the value of /3 using algebraic arguments along with some 
heuristic hypotheses [TO]. At this point, it is important to notice that it is also possible to compute 
/3 using the simple limit process. This was already noticed by Gurarie and Ludwig in the context 
of the replica approach [21]. Using the parametrization a; = 2 + e, we find 

c/2 5 

/^percolation = ~ 1™ 7 ;^=~o' (1.13a) 

€— >0 H-i 5 — Z o 

c/2 5 

^polymers = -1^7 ^ c ' (1.13b) 

which indeed are the expected values [21 HO]- In this sense, it is possible to understand the structure 
of a logarithmic CFT as a limit of non-logarithmic CFTs. We suspect that the values of indecom- 
posability parameters can be inferred in a similar fashion for any LCFT. To see this, we turn to a 
slightly more complicated example. 



1.3 Generalization to other LCFTs and general formula for 13 

Let us consider a generalization of this argument to a more complicated case. We focus in this 
paragraph on the theory of symplectic fermions [271 [2H| that describes the scaling limit of dense 
polymers on the lattice. This theory has c — —2 (x = 1), and is particularly simple as there 
is no interaction. The Jordan cell structure can be understood in terms of free fermions, and the 
indecomposability parameters can be, in principle, computed using this free fermion representation. 
This theory was also considered in order to test the Nahm-Gaberdiel-Kausch algorithm to compute 
fusion products between indecomposable Virasoro modules [12]. 

We choose here to focus on a concrete example: this theory is known to have a Jordan cell at 
level 3 with equations 

Lo(/. = 30 (1.14) 
LqiP — Sifj + (j) 

A^L\-2L^2, 

where the parameter /? = {4>\''l') = — 18 is the logarithmic coupling associated with this Jordan cell 
[12]. Note that we chose a different normalization convention for the operator A than eq. ()1.4|) 
in order to match [12]. When we think of this theory as the continuum limit of a XX spin chain 
with quantum Uq=i{5\.2) symmetry with an even number of sites [15j . the conformal dimensions 
that appear in the spectrum are hi^i+2j with j € N. We shall come back to the precise nature 
of the scaling limit of such lattice models in the next section. All we need for what follows is to 
identify the conformal weights of the fields ^ and -if) in the spectrum. We find — /11.5 = 1 and 
ft-i/i = /ii.7 = 3. It is important to identify the operators in spectrum as we are interested in small 
perturbations around c = —2. Let us now consider a conformal field theory slightly deformed from 
c = — 2, with X = 1 + e. Within this generic (non-logarithmic) CFT, the OPE of a generic operator 
$^ with itself reads 



<i>.(.)<i>.(0) ^ 



m + \zdm + a(-2)22^_2^(o) + a(-i^-i)z2L^ie(0) 



(1.15) 
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where the a coefficients are fixed by conformal invariance and are diverging as e — ^ 0: a*- — 
+ i±|^ + 0{€) and a(-i'-i) = + + ^(^)- ^ote that we showed only the channel of ^ 
(with dimension = /ii 5 = 1 + |f -l '''^^^2^ 
= (L2_^ _ 2L_2)e 



1 + |e + Oie^)) on the right-hand side. Let us now introduce the field 



^2h — hc 



(1.16) 

Remark that we got rid of one of the diverging terms this way as Orcg + = 

+0(6). At this point, we have no choice but to admit that there exists another field in the 
theory with dimension 3 at c = —2 to cancel the last diverging term. As we already discussed, this 
field has to be $1^7 with conformal weight = hi^r = 3 + 3e + O(e^). The OPE thus reads 

+a,egz''«+iL_2e(0) + z''*'-i$,_7(0) + ...]. (1.17) 

We define the field ipiz) as 

v(-l:-l) 



*1,7(^) 



Pie) 



Tpiz) — a 



(-1,-1). 



where /3 
)ecomes 

^h{z)M0) 



(010) 



/i^_ftj_2"- We finally notice that a' ^\h^ — — 
becomes regular at c = — 2 and involves logarithms 



{zl (1.18) 
= -^+ 0(e), so that the OPE 



^2/1-1 



m + ^om + ^^z^L^2m + ^z\m + m log^) + 



(1.19) 



One can check that the operators t/i and (j) defined this way satisfy the usual OPEs for logarithmic 
operators. In particular, it is straightforward to show that 



{<i^{z)m) 



/3 



1,7 



(1.20) 



where 



Bn — lim /3(e) = — lim - — 



m) 



= -18, 



(1.21) 



as (010) = 27e + C'(e^). We find the same /3 parameter as [12] but with a different (technically sim- 
pler, but less rigorous) argument which only involves computation of a few Virasoro commutators. 



We now turn to a more general LCFT with central charge c — 1 



Using the previ- 



xo(xa + l) ■ 

ous results, we conjecture that the indecomposability parameter (5 for a generic Jordan cell with 
structure (|1.3p , can be computed from small deformations around this theory x — xq + e as 



B- 



lim 



e=0 



h,! 



[h^p - h^) 



E=0 



(1.22) 



where n — {h^ — h^)\^^Q- We will show in the two next sections that eq. (|1.22p is consistent with 
the previous studies and agrees very well with numerical results. The problem now reduces to 
identifying Jordan cells of a given theory. This problem is fairly well understood, and we shall now 
turn to concrete lattice examples to illustrate this. 
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Figure 1: Graphical representation of the Temper ley-Lieb-based vertex model. The lattice consists 
of alternating arrows going up for i even and down for i odd, where z = 1, . . . , L = 2N corresponds 
to the horizontal (space) coordinate. The system has free boundary conditions in the horizontal 
direction and periodic in the vertical (imaginary time) direction. We choose each vertex according to 
its probability; this draws a dense loop configuration on the lattice. Each closed loop carries a weight 
n ~ q + q^^ ■ In the supersymmetric language, the alternating □, □ representations correspond to 
a lattice orientation, conserved along each loop. The system is isotropic when pA = Pb, while the 
transition occurs when pA = 1 — Pb ■ 

2 Lattice models and algebraic considerations 

Although the systematic study of Virasoro indecomposable modules can be analyzed on the Virasoro 
side in an rather abstract way [121 114j , it is also instructive to analyze how indecomposability arises 
directly from lattice regularizations [T21 [IS] ■ The structure of the Virasoro staggered modules can 
be predicted from the analysis of the projective modules of associative lattice algebras such as 
the Temperley-Lieb algebra [13 [Ul HH] . These results are of course consistent with those of the 
Virasoro-based approach. This section follows closely the results of Ref. [15] . 

2.1 'Dense' LCFTs from Temperley-Lieb algebra, XXZ spin-1/2 chain, 
dense loops and supersymmetry 

2.1.1 Temperley-Lieb algebra 

Underlying most of the models we shall consider is the Temperley-Lieb algebra TL2N{q)- The 
algebra TL2N{q) defined on L = 2N {N G N/2) strands consists of all the words written with the 
2N — 1 generators (1 < « < 2N — 1), subject to the relations 




o(N-j|>2) 



(2.1a) 
(2.1b) 
(2.1c) 
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with 

n^q + q-\ (2.2) 

We also define q = e'^ and 7 = ^j^. The action of as an operator acting on 2N strands can 
be naturaUy represented graphically. This constructs what we shall refer to as loop or adjoint 
representation. The standard modules of the TL algebra are well known, and have dimensions 

where in the geometrical language, 2j is the number of 'through lines' or 'strings' that propagate 
along the imaginary time direction. We have j G N, restricted to the condition + j G N, so if 
is a half-integer (L = 2N odd), so must be j. These standard modules are irreducible for q generic 
(i.e. not a root of unity). 

We consider a two-dimensional model defined by the transfer matrix 

N-l N 

T=l[{PB + {l- Pb) 62.) n ((1 - pa) + PA 62.-1), (2.4) 

i=l 1=1 

which acts on a given TL module. This definition is valid for L = 2N even but it can be readily 
adapted to an odd number of sites. We will mainly work with three different representations: 
geometrical (adjoint), 6- vertex, and supersymmetric. Using the geometrical (loop) representation 
of TL2n{ci), we obtain a dense loop model, where each closed loop carries a weight n (fugacity). 
A graphical representation of this vertex model is given Fig. [1] It sometimes helps to think in 
terms of a Q-state Potts model whose high-temperature expansion consists of drawing dense loop 
configurations with fugacity n = \fQ. 

In the strong anisotropy limit ^ with pa/(1 ^ Pb) fixed, we can extract the Hamiltonian 
of the equivalent one-dimensional quantum system. It reads, up to an irrelevant constant, 

N-l N 

iJ=-£^ 62. -£-1^62.-1, (2.5) 



where e — \J PaI (1 ^ Pb)- The system is isotropic when "pA ^ Pb, while the transition occurs when 
PA = 1 ^Pb- Hereafter, we will always consider the case e — I. 

Finally, we recall some useful relations for this integrable model. Using Bethe ansatz, one can 
show [53] that the mean value of the TL generators on the groundstate is 

2 /■+°° da; 1 ^2 g-j 

°° J-oc coshTra; cosh 27a; — 2 cos 7 

Meanwhile, the dispersion relation of the quasi-particles is e ~ fc^^j, where the sound velocity 

reads 

TT sin 7 

VF = 1. (2.7) 

7 
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2.1.2 6- vertex or XXZ spin chain representation 

Another natural representation of the Temperley-Lieb algebra is provided by the 6-vertex model. 
We write n = e^'' +e~^'^ and q = e^^ . The Hamiltonian limit of the 6-vertex model is the XXZ chain, 
with Hilbert space "Hxxz = (C^)®^^. We will focus on this limit hereafter. The Temperley-Lieb 
generators in this representation read 



/ 







q-^ 


-1 





-1 


q 





\Q 





/ 



(2., 



where we have used the basis {|tt) , Iti) , lit) , lit)} of (C^)®^, the tensor product of the Hilbert 
spaces of the sites i and i -I- 1. We can check that the generators indeed satisfy the TL algebra, 



with n = q 
operators 



We can express the Hamiltonian of the XXZ chain in terms of the usual Pauli 



2N-1 

2 ^ 

i=l 



q + q" 



'i+i 



'i+l 



K-Cr|^)-iVcOS7. 



(2.9) 



Note that this XXZ Hamiltonian (or equivalently the transfer matrix of the 6- Vertex model) com- 
mutes [30] with the generators of the quantum group Uq{5i2)- 



2.1.3 Supersymmetric representation (SUSY) 

We introduce in this paragraph a supersymmetric formulation of our model, which provides another 
natural representation of TL2n{<i) [31] • We consider that each edge of our two-dimensional lattice 
carries a Z2 graded vector space of dimension n + m\m. We choose these vector spaces to be the 
fundamental □ of the Lie superalgebra gl(n + m\m) for i odd (corresponding to down arrows of 
Fig. HI), and the dual □ for i even (up arrows). The transfer matrix (or the Hamiltonian) then acts 
on the graded tensor product V. — {D D)^^ . The TL generators are simply quadratic Casimir 
invariants, providing a natural generalization of the Heisenberg chain to the Ql{n + m\m) algebra. 
We can check that a loop expansion of the transfer matrix yields a dense loop model with a weight 
stT I — n + m m — n for each closed loop as expected. There is a continuum quantum field 
theory associated with this spin chain, which turns out [311 to be a non-linear cr-model on complex 



projective space 



^n-t-m — l|m 



U(TO-t-n|TO)/(U(l) X U(m- 
2.1.4 Continuum limit: the generic case 



l|m)) at topological angle = tt. 



Let us also say a few words about the continuum limit. The XXZ chain with q = e*''', j — 7t/(x+1) 
is described by a CFT [30j with central charge 



0= 1 



6 



x{x + 1) ■ 



(2.10) 



Consider the Verma modules V^.s, spanned by the action of the Virasoro generators Ln with n < 
on the highest weight state with conformal dimension hr^s given by the Kac formula 



[{x + l)r - xsf - 1 

4:X{X+1) 



(2.11) 
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When we take the scaHng limit of the XXZ chain, the conformal weights that occur in the spectrum 
are hi^i+2j in the Kac table, where j — n, n S N if L even and j = n/2, n £ 2N + 1 if L odd. For q 
generic (not a root of unity), there is a single null vector in the Verma module Vi,i+2j at conformal 
weight hi^-i-2j that must be set to zero in order to obtain a simple (irreducible) module. Hence, 
we define the standard (Kac) modules rj — ri^i+2j — Vi,i+2j/Vi.-i_2j which are irreducible for q 
generic. The charactejfl of the module rj reads 

«'»l,l + 2j-c/24 _ /ti i_2j-c/24 

K, = Tr, ,--^/- = 1 . (2.12) 

where P{q) is the inverse of the Euler partition function, and is related to the Dedekind r] function 

oo 

P{q)^l[i^-qn = q-'^^'viq)- (2.13) 

n=l 

The partition function of the sector of spin reads 

oo 

^XXZ ^ qL„-c/2i ^ j^^^ (2.14) 

j=\s.\ 

the global partition function of the Uq{sl2) XXZ spin chain is then readily obtained 

+ 00 oo 

Z^^^^ Y: Zf^-^Y.(^j + l)K,. (2.15) 

S^ = -oo j=0 

The partition function of the superspin chain with gl(n + m\n) symmetry is given by a similar 
expression where one replaces 2j + 1 by the dimensions of the irreducible representations of the 
commutant An+m\ni'^N) of TL2n{q) in the SUSY representation I32j. The algebra An+m\n is in 
fact much larger than g[(n + m\n). 

These results were obtained for q generic but are supposed to remain correct even when g is a 
root of unity even though the Virasoro standard modules are no longer irreducible in general. 



2.2 'Dilute' LCFTs from the integrable 0{n) model 

We will also consider a fundamentally different version of the previous LCFTs using the integrable 
dilute 0{n) on the square lattice. We shall refer to these theories as 'dilute' as opposed to the 'dense' 
ones based on the Temperley-Lieb algebra. This denomination obviously refers to the dense or dilute 
nature of the underlying loop gas. Therefore, we describe in this section the 0{n) model defined on 
an annulus of width 2N. It corresponds to a dilute loop model where closed loops carry a weight 
n; we shall focus here only on the dilute phase. This model also possesses a dense phase which is in 
the same universality class as the dense loop model. Note that the case n H> is obviously relevant 
for the physics of polymers. In terms of spin chains, it is described by a 5* = 1 C/g(s[2)-invariant 
chain where the states Sz — are viewed as occupied by parts of loops and 5*2 = as empty. 

^The modular parameter q in characters and partition functions has of course nothing to do with q = c^'' that 
parametrizes the weight of closed loops in the Temperley-Lieb algebra. We keep the same notations and hope this 
will not confuse the reader. 
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P3 



P7 



Figure 2: Plaquettes and integrable weights for the dilute 0{n) model on the square lattice 

This model also corresponds to osp(n + 2m|2m) (super)spin chains and to non-linear sigma models 
with supersphere target space 52m+n-i|2m ^ 0Sp(2to + n|2m)/OSp(2TO + n- l\2m) [2T]- There 
is a dilute version of the Temperley-Lieb algebra behind all these models. We will not go into the 
details of these different formulations here but only describe the geometrical setup that will allow 
us to measure indecomposability parameters in the next section. 

2.2.1 Lattice model 

Our starting point is the expression of the integrable version of the dilute 0{n) model defined on 
the square lattice. The R matrix is the sum of 9 plaquettes represented graphically in Figure O 



(2.16) 



where j G {1, . . . , 2N — 1} labels the sites. It satisfies the Yang-Baxter equations for the corres- 
ponding integrable weights [35] 



Pi 

P2 

Pi 

P6 
P8 
P9 



1 



sinwsin(3A — u) 
sin 2A sin 3A 
sin(3A — u) 



sin 3A 
sin It 



sin 3A 

sinusin(3A — u) 



P3 = 
P5 = 

sin 2A sin 3A 
sin(2A — u) sin(3A — u) 
sin 2A sin 3A 
sin u sin(A — u) 
sin 2A sin 3A 



(2.17a) 
(2.17b) 
(2.17c) 
(2.17d) 
(2.17e) 
(2.17f) 
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where n — 



— 2cos4A is the weight given to every closed loop. Note that 



R,(m = 0) = Ol + Ol + Ol + 0\ = Id. (2.18) 

One can extract the corresponding ID Hamiltonian H = — ^ using an expansion in the 
spectral parameter u. We find that the interaction between the sites i and « + 1 reads 

= -^0\ - cotan3A(02 + o\) + -KrzkOl + 0\) 
sm zA sm 6\ 

+ -^71(0.' + - (cotan2A + cotan3A)Of - . "'"^ ^. O^ (2.19) 
sin2A sin2Asin3A 

In order to obtain an integrable Hamiltonian in the case of open (free) boundary conditions, we 
would need to consider K matrices at the edges [31]. Instead, we here consider the simpler (non- 
integrable) Hamiltonian R = — J2 ^i-i+i- what follows, every reference to the 0{n) model 
corresponds to this Hamiltonian. The Fermi velocity and the mean value of on the ground 

state are |35j 

VF = ^, (2.20) 

and 

tanhSAfc cosh(5fcA — /ctt) cosh fcA 

d/c ^— 1 . (2.21) 

smh KTT 



2.2.2 Continuum limit 

The CFT describing the scaling limit of dilute loops has central charge [3S] 

c=l-3 — — , (2.22) 

ZATT 

with n = — 2cos4A and A G [7r/4,7r/2]. The conformal weights appearing in the spectrum lie now 
in the first row of the Kac table /ii+2j,i- The trace over the standard module with 2j through lines 
reads 

K, ^ Tvr^ = i . (2.23) 

As in the dense case, the partition function of the spin-1 XXZ chain is just the sum of these 
characters with a degeneracy 2j + 1. In the case of osp(n + 2m|2m)-invariant superspin chains, the 
degeneracy 2j + 1 is replaced by the dimension of the corresponding irreducible representation of 
the commutant 'B„+2m|2m(2A^) of the dilute Temperley-Lieb algebra |i32j . 



2.3 Indecomposability and lattice Jordan cells at q root of unity on a 
concrete example: the logarithmic Ising model £7W(3,4) 

We now turn to the study of indecomposability when g is a root of unity. In this case, the standard 
modules of the Temperley-Lieb algebra still exist but may no longer be irreducible. Read and Saleur 
[TS] showed in details that the structure of the projective modules of the Temperley-Lieb algebra 
parallels that of several staggered Virasoro modules in the continuum limit. The general algebraic 
structure of the XXZ spin chain was illustrated in terms of 'staircase diagrams' as a representation 
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(bimodule) of TL2n{<i) ® [^^(sb), or TL2N{q) ® An+m\n{'^N) in the SUSY case. In this paper, we 
win pass directly to the continuum theory and only mention the important algebraic lattice results 
when necessary. We refer the interested reader to Ref. [TS] for more details. 

Let us consider the case q = which corresponds to a dense CFT with central charge c=\. 

Loop expansion of the partition function yields a loop model with fugacity n = this is of course 
the celebrated (logarithmic) Ising model £A1(3,4) [TB]. The Kac formula at a; = 3 reads 

Ks - ['^'^f '\ (2.24) 
and the values appearing in the spectrum of the t/g (5(2) -invariant XXZ spin chain at g = e^'^l'^ are 

fti,i+2i = • (2-25) 

We will focus the continuum limit of a spin chain with an even number of sites so that j G N; the 
case 2 = n/2, n g 2N + 1 is treated in a similar way. Note that we have hr^s = h-r^-s = ^t'+3,s+4- 
The character of the standard module rj reads 

K, - Tr,^ = 1 ^ , (2.26) 

whereas the characters (37) of the simple (irreducible) Virasoro module are 

^(24n-l+24p)^/48 „ j^(24n+17+24p)^/48 

Xj=4p = ^7^^ ' (2.27a) 

ni[-2p.-l] '^^> 

(24n+5+24p)^/48 _ (24n+ll+24p)^/48 

X.=4p+i - E 17^? ' (2-27b) 

(24n+ll+24p)^/48 „ (24n+29+24p)^/48 

X,=4p+2 = E ^- , (2.27c) 

ri^[-2p-l,-l] 

(24ri+17+24p)^/48 _ (24ri+23+24p)^ /48 

X,=4p+3 = E ^7^^ ■ (2-27d) 

n^[-2p-l,-l] 

We can see using these character identities that the standard modules rj, j e N are no longer 
irreducible. The decomposition onto the simple characters indeed yields 

Kj=ip = X4p + X4p+3, (2.28a) 
i^j=4p+i = X4p+i+X4p+2, (2.28b) 

ifj=4p+2 = X4p+2+X4p+5, (2.28c) 
ifj=4p+3 = X4p+3 + X4p+4- (2.28d) 

Let Rj be the simple Virasoro module with conformal weight /ii^i+2j and character Xj- From lattice 
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algebraic considerations, we know that the standard modules rj must thus have the structure 



R 



•i+3 



R, 



j = (mod4) or j = 2 (mod4) 
j = 1 (mod4) or j = 3 (mod4) 



R 



The arrows represent again the action of the Virasoro algebra. The bottom submodule is invariant 
while it is possible to go from the top to the bottom acting by some element of Virasoro. The top 
simple quotient must have a smaller j number than the bottom |15j . 

Using the standard modules as elementary bricks, it becomes quite easy to construct the 
staggered modules of the theory using the knowledge, again from the lattice, that they must be 
diamond-shaped. For instance for j = 2, we expect 



R2 

\ / 
R2 



(2.29) 



with character P2 = 2x2 + Xi + X5- This is the only gluing of standard modules that respects the 
fact that the conformal weights must be increasing to the right. For a given theory, there is always 
a unique way to construct a given diamond-shaped module in terms of simple modulef[^ using only 
the structure the standard modules. 

Using this method, it becomes quite straightforward to guess the general structure of the Vira- 
soro staggered modules appearing in the theory 



Ro 



Ri 



R3 



R2 



j = 0, 



Ri 



\ 



R 



R 



J+3 



j = (mod 4) and j > 0, 01 j = 2 (mod 4), 



(2.30) 



J-3 



Ri 



Ri 



Ri 



R 



j — 1 (mod 4) and j > 1, or j = 3 (mod 4). 



^To be more precise, we mean here that for a given module Pj of a given theory (dense or dilute), the values of 
the subscripts j of the four simples Rj in the diamond Pj are uniquely fixed. 
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Of course, it is also possible using Temperley-Lieb representation theory to check that a similar 
pattern arises from the lattice. Note also that it may happen that the standard module rj remains 
irreducible for q a root of unity. For example, in the Ising model with L odd, one finds that 
ifj=4p+3/2 = X4p+3/2 so that what we shall call staggered module in this case is just the simple 

module itself Pj=4p+3/2 = rj=4p+3/2 = Rj=4p+3/2- 

Let us summarize what we have learned so far concerning the structure of the scaling limit of the 
XXZ spin chain at q = e^'^/^ with an even number of sites. All the fields of the theory are organized 
into staggered modules given by (|2.3U|) . j £ N. Recall that each simple module Rj corresponds to 
a field with conformal weight hi i^2j- There is a Jordan cell in Lq for every staggered module with 
a diamond shape. For such generic staggered module, we note the basis fields 

Rj V^^'^ 
^ \ / \ 

V,= R,, R,, = , (2.31) 

R, 

with ji < j < 32 so that /ii,i+2ji < /ii,i+2j < /ii,i+2j2- We define the logarithmic coupling /3i,i+2j 
of this Jordan cell as 

/?i4+2, = (0^^'V^^'^)- (2.32) 

Note that we use the same Kac labels for j3 as those of the conformal weig ht of i/i^J'). 

We are now ready to compute this indecomposability parameter using eq. (jl.22p along with 
— ft.i,i+2ji and — /ii^i+2j2- Let us again illustrate this on a concrete example. The first 

Jordan cell occurring in the spectrum has j = 2 and its structure is described by eq. ()2.29p . Defining 

a; = 3 + e, it is straightforward to compute hi^^ = 5 + f + ^i^^) = 1) ^"^^ ^i,5 = I + T + '^('^^) 

(j = 2). At e = 0, the Lq operator expressed in the basis {(j)^'^\ip'-^'>) reads 

Lo - ( ) . (2.33) 

We remark that we can readily find the relation between 0^^^ and (up to an irrelevant global 
normalization factor) because we know that (j)^^^ must be a null-vector 

c^('H^) = iL-2-lLU)e'\z). (2.34) 

In general, the operator A is either known from general formulas or computed numerically using 
the Virasoro algebra. Straightforward commutations of the X„'s modes yield the final result 

£^0 hi.5 - hi.3 - 2 24 ^ ^ 



2.4 General structure of the staggered Virasoro modules 

It should be clear that the path followed in the previous subsection can be extended to all our 
models. In this paper, we focus on the dense and dilute versions of the minimal logarithmic models 
£A4{l,p) and £M{p,p + 1), p & N. The former choice corresponds to a; = 1/p and the latter to 
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X — p, with q — e"/*^^"*"^^. The structure of the standard modules in all cases can be inferred from 
characters identities, and the staggered modules can then be built using the standards as elementary 
bricks. Using character identities, it is possible to convince oneself that the two following statements 
should hold 



• the CA4{l,p) theories (x 
with X = p, 



1/p) have the same structure as their 'dual' £M{p,p+l) theories 



• a 'dilute' LCFT based on the 0{n) model has the same structure as the 'dense' LCFT with 
the same loop fugacity n = q + . 

We say that two theories have the same structure when the Vj modules in both theories have the 
same expression in terms of the simple modules Rj, that is to say that they are characterized by 
the same values of ji and j2 in eq. (|2.3ip . This docs not mean that the modules are the same, in 
particular, they may be characterized by different indecomposability parameters, and the simple 
modules Rj are completely different as they are over different algebras. For example, the staggered 
modules in the theory with x = 1/3 have the same structure in terms of simple modules as those 
of the Ising model, given by eq. (|2.30p . 

Therefore, everything boils down to the study of the staggered module^ of C^A {p,p+l) theories 
with x — p. For such theories, we find that the staggered modules have the following form 



-Pi 



R 



R 



j+p 



(mod and j > 0, 



Rj 



\ 



Rj— 



a" -2 



R 



j =\ (mod 2±i) and j > 1, 



Ri 



Rj 



\ 



Rj-f 



R-j 



Rj 



R,+i j = ^ (mod £±i) and j > p - 1, 



J^i (mod£±i). 



Once again, each diamond module corresponds to a Jordan cell for Lq involving the null-field (t>j{z), 
with logarithmic coupling /3i,i+2j {resp. /3i+2j,i) in the dense {resp. dilute) case given by eq. (|1.22p . 
For j < p — 1, only the top standard module of the diamond remains in the staggered. 

^Note that what we call staggered modules here are nothing but the 'scaling limit' of the Temperley-Lieb projective 
modules that arise in lattice models. They may or may not be indecomposable diamonds depending on the value of 
the 'spin' j. 
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3 Numerical measure of indecomposability parameters from 
lattice models 



While the analysis of symmetries of the lattice models provides information about the general struc- 
ture of the Virasoro indecomposable modules, getting more detailed information about the action 
of the Virasoro generators in these modules — such as the numerical values of the indecomposability 
parameters — is more challenging. 

There are many difRculties to overcome in this kind of analysis. One of the most important 
ones concerns the proper normalization of the state 0, which obeys {4>\4>) = 0. For c = and the 
Jordan cell of the energy-momentum tensor, Ref. 20] used a trick that led to the determination of 
b — —5/8 for percolation and h — 5/6 for self-avoiding random walks (dilute polymers). We follow 
here a different route, and propose a more general method which allows us to study other Jordan 
cells with an A operator more complicated than L_2. The method turns out te be quite accurate so 
that indecomposability parameters can be determined numerically, with almost the same precision 
as for critical exponents. 

3.1 Jordan cells and lattice scalar product 
3.1.1 Spectrum and Jordan cells 

We consider a generic Jordan cell at level h in Lq. As before, we normalize our states such that in 
the basis {4>,tp — A£^), Lq reads 

io^(j I). (3.1) 

It is well known that the operator Lq can be related to the scaling limit of a Hamiltonian at a 
critical point. For a system of length L — 2N, we have 

H = Es + E^L -f ^ (Lo - ^) + 0{L-\ (3.2) 

where Es is a (non-universal) surface energy and E^ is the bulk energy density. Hence, we 
expect non-diagonalizability also in our lattice Hamiltonians (or transfer matrices) that mimic the 
continuum limit behavior. The central charge and the critical exponents of a given model can be 
readily measured numerically using finite-size corrections \n L — 2N to the eigenvalues of H. If we 
note Eq the energy of the fundamental and E^ that of a given excitation, we have the well-known 
relations [35| 

En{L) = Eq{^)L + E^ -^^ + ... (3.3a) 
E4L)-Eq{L) = E^ -E^ + ^h^ + ... (3.3b) 

with most of the times E^ — Eq . We here assumed that all the conformal dimensions were 
positive; if that is not the case, one introduces the usual concept of effective central charge Ccs — 
c — 24/imin. Using these formulae, it is then a simple matter to identify the different eigenstates of 
the Hamiltonian which corresponds to {ip), \(f)) or |^) in the continuum limit. To do so, we use the 
Arnoldi algorithm [39^ to get the eigenvalues and the corresponding Schur vectors for the lowest 
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excitations of the spectrum. We then apply some variant of the Gauss- Jordan algorithm to put the 
reduced Schur (upper triangular) matrix into Jordan canonical fornH. 

Let us suppose that we identify a Jordan cell in the Hamiltonian which corresponds to an energy 
E{L). We normalize the states to prepare the comparison with CFT. In the basis {|'/>^^'') , 
the Hamiltonian for a system with L = 2N sites reads 

where vp is the Fermi velocity and h^^^ = ■:^^{E{L) — Eq{L)). Note that we have limi_j.oo /i*-^^ = h, 
so the matrix expression in eq. (j3.4p goes to Lq in the continuum limit. 

3.1.2 Lattice scalar products 

In order to measure P = {iplcj)), we first need to define a 'scalar product' that goes to the Virasoro 
bilinear form in the scaling limit. The construction of lattice 'scalar products' which go to the 
Virasoro form in the continuum has already been studied in great details in Ref. ^U\. All scalar 
products must be chosen such that = Lq in the underlying CFT. This means that we want the 
Hamiltonian H to be hermitian for these scalar products. Of course, as we deal with non-unitary 
theories, there may be negative-norm states so what we call here scalar product is nothing but a 
sesquilinear form. Obviously, the non-hermitianity of H for the usual scalar product was the reason 
for its non-diagonalizability in the first place. Let us recall the expression of the scalar product for 
the different representations of the (dense/dilute) Temperley-Lieb algebra: 

XXZ: The scalar product is the usual bilinear form on C without complex conjugation, that is, 
treating g as a formal parameter. For example, on L ~ 4 sites, the vector |0) — |ttii) + Q Itttt) 
has norm {(j)\cj)) = 1 + q^. Note that if we had considered the usual scalar product on C, we would 
have found 1 + \q\'^ instead. 

LOOP: The correct scalar product is obtained gluing the mirror image of the first state on top of 
the second one. Each closed loop carries a weight n = q + q~^. This is of course the usual form used 
in Temperley-Lieb representation theory. For instance, the scalar product between the two states 
|a) = Ik^) and |/?) = |^ is (a|/3) = = (/\) = n. The case with a non-zero number of 

strings 2j is treated in a similar fashion. Finally, in the case of the dilute 0{n) model, the scalar 
product between two states is chosen to be zero if the empty sites (marked as dots) are not the same. 
For example for i = 6, * |^-' * u •) = (/^^^ _ |vy • u •) = whereas * * I*-" * ^ *) = 0. 

SUSY: We use the usual scalar product in Fock space. There are negative norm states because of 
the use of the dual representation. For example, let us consider the s[(2|l) case still on L = 4 sites. 
A precise definition of this chain will be given in the following. The important point here is that 
each site must be occupied by one particle which can be either a fermion = {—iy~^^Sij, 

or a Schwinger boson [6^,^,^]^'] — ^ij^aa', with a e {tii}- Let us consider the state |0) = 
b\^f2bl^b\^\0). Its norm is = (0| 64^63^/26it&i^/2&34,&4-f- 10) — —1 because of the fermionic 

operator of the dual representation which satisfies {/2, = — 1- 

*This step involves choosing a rule to declare that two numbers Ai and A2 are equal if and only if |Ai — A2I < £ 
with typically e ~ 10~*. 
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3.2 Virasoro algebra regularization on the lattice 

There is a last difficulty that one must tackle in order to define a proper version of j3 on the 
lattice. Remark that the Jordan cell in (|3.4p is invariant under a global rescaling of the basis 
states \(l>'"^'^) — > al^^^-*) and jV'''"'^^) a\'ip''^^). Unfortunately, the scalar product between the 
states is not invariant under such transformation ^f/j^^^ — > |ap Hence, we need 

to normalize the state l^*-^-*) so that it goes precisely to |0) = A \£) in the continuum limit. 

Let |(/''"^-') ~ a If^'"^-*) and jV''-"^'') = a where goes to |0) when L — > oo. If we knew 

a, we would be able to compute — 7> /3 as L ^ oo. Note that is a null state so that 
we cannot simply normalize it, we thus need to find another way to get rid of this global factor 
a. This is achieved using a regularization of the Virasoro generators on the lattice. Following Ref. 
PT| . for a general critical Hamiltonian 

L-l 

= (3.5) 
with Fermi velocity VF^ we define a lattice version of the L„'s 



TT 



1 /, , / niT:\ 1 , , . f nin 



/ (/li - /loo) cos — — + ^ > [/ii, ft-i+i] sin 
i—i ^ ^ ^ i—i 



Sn,0, (3.6) 
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where /loo is the ground state expectation value of hi. It is possible to show that such expression 
provides a good way to define the Virasoro algebra on the lattice in the case where hi = is 
a Temperley-Lieb generator. In particular, it is possible to measure the central charge or scalar 
products of the continuum limit through the computation of Virasoro commutators. We shall 
admit that this equation remains correct even if hi is a generator of the dilute Temperley-Lieb 
algebra. Note that these generators do not exactly satisfy the Virasoro algebra in the continuum 
limit because of anomalies; the commutators of the scaling limit do not coincide in general with 
the scaling limit of the commutators due to extra couplings to 'non-scaling states' [H]. While the 
problem can be solved using a double limit procedure, in practice, the anomalies induce extremely 
small errors to the lattice measurements. For practical purposes, formula p.6p can thus be used 
naively, even when commutators of multiple actions of Virasoro generators are involved. 

Using formula p. 61 we are thus able to construct a lattice version A^-'^^ of the A operator that 
links (f> and ^. Moreover, the state is readily identified in the spectrum. We normalize it such 

that (C'^^^IC*'^^) = 1- If we assume that \4>''^^) = A'-^^ is a correct lattice version of \(f>), we 

are now ready to compute /3. Gathering all the pieces, we define a lattice version of /? which does 
not depend on a 



(3.7) 



When L — !> oo, we check that 



Up |/J,(L)M(L)t(L)\ 1^ . , 

hm = lim ' ' ^'LjLaLx = (^^^^I'^^'^O = ^^l'^) = (3.8) 

L->oo L^oo \a\ {'Ip'- '\(p'- >) L^oo \ I 

We summarize our method to measure /3 by the following steps: 
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1. Using exact diagonalization methods, find a Jordan basis for the first few excitations of H on 
L^2N sites. 

2. Identify a Jordan cell in the spectrum of H and normalize the states like in eq. p.4p . 

3. Also identify the state and normalize it such that = 1 for the lattice scalar 
product. 

4. Using Virasoro generators on the lattice (|3.6p . construct the operator A^^\ 

5. Compute ^(^^ using eq. (|X7|) . 

The value of the indecomposability parameter /3 = MraL^aa /S'^-* is then computed using an extra- 
polation Z?'^) = (3 + A/ L + B / L'^ + . . . We find numerically that does not depend on the chosen 
Temperley-Lieb representation (Loop, XXZ, or SUSY). However, it does depend on how we realize 
the operator A on the lattic^, this is why we were able to improve the results of Ref. [^U] . 



4 Numerical results 

We present in this section our numerical results for indecomposability parameters in different ex- 
amples of LCFT. 



4.1 The case x = 1 (c = —2): Dense Polymers 

We begin with the case of the XX(Z) chain of even length 2N, with q — i so c = —2. The expansion 
of the partition function in terms of dense loops describes dense polymers as the weight for a closed 
loop is n = 0. This theory is relevant for the description of spanning trees [12] (up to a duality 
transformation). It is also related with the description of abelian sandpile models, although in the 
latter case different values of the indecomposability parameters have been found [101 HI] ■ In the 
SUSY language, this model also corresponds to g[(m|m) (super)spin chains and to non-linear sigma 
models with target space 

^pm-il™ at = TT pri. For m = 1, we get the 0[(1|1) spin chain which is 
a free fcrmion system. Indeed, at m = 1 the Uq{sl2) XX spin chain and the supersymmetric 0[(1|1)- 
invariant chain coincide and everything can be reformulated in terms of free-fermion generators 
obeying {fi, fj} — and = {—ly^^Sij. Within this representation, the Temperley-Lieb 

generators read 

e. = (// + ./^i)(/. + /.+!)■ (4.1) 
The corresponding continuum limit is a symplectic fermions theory [271 128] with action 

^ = ^ y d^eO'^t (4.2) 

This theory is probably one of the best understood LCFT. There are 4 different fields at level which 
are organized into a diamond indecomposable module of g[(l|l). All the Virasoro staggered modules 

^In general, there are infinitely many ways to realize the operator A on the lattice. In the case of yl = -L_2 acting 
on the vacuum for example, one can use the trousers trick [20], or the Koo-Saleur formula 1 13.61 1 directly A = L^^', 

or even A = i*!^' -I- aL^^^ where o £ C. The values of computed using these different lattice realizations are 
different, although we expect them to yield the same result in the limit N — >■ oo. 



22 



can be constructed as fermionic excitations of these states at level 0. A few indecomposability 
parameters have been computed using the Kausch-Gaberdiel algorithm [T^]. One can also readily 
compute the same parameters using the free fermion representation. 

We show here how to implement our lattice approach to this model. The Kac formula at a: = 1 
reads 

(2r-sf-l 



-r,s- 3 , (4.3) 

and extending the generic results to the case q = i, we find that the partition function of the 
spin chain is 



Z = J2i2j + = ^j(2x,M +X,-i,i)- (4.4) 

j=0 j=l 

Recall that Xj,i — Xi.i+2j is the irreducible character of the Virasoro simple module with h = 
hj^i = /ii,i+2j- Using the general results presented in section [2.41 with p = 1 and j = ^(mod 1), 
we see that there is no indecomposability for L odd {j half-integer) so that the lattice Hamiltonian 
remains fully diagonalizable in this case. For j e N, the staggered Virasoro modules have the 
following subquotient structure and a basis in terms of fields 

-Pj^ Rj-i Rj+i - . (4.5) 

\ / \ v/ 

The modules Vj are completely characterized by the logarithmic couplings /3i^i+2j, and each module 
corresponds to a Jordan cell 

Lo0(^') = /ii,i+2,(/'(-'') (4.6) 

Note that we choose a different convention for the operator A than eq. (|1.4p in order to match [T^ 
as we normalize it such that Aj = L^Si + . ■ ■ Using this convention and eq. (|1.22p , we were able to 
conjecture a general formula for the indecomposability parameter 

/3M+2,=-lim^ f'^^^'^^ =- [^^^~_f^' (j-l), jeN\{0,l}. (4.7) 

/ll,l+2j - /il,l+2(j-l) - (j - 1) 4J ^ 

We checked this equation up to j = 16 using eq. p.22p . Note also that in other cases can similar 
explicit formulae be obtained from eq. (|1.22p ; we shall report on this point and compare our results 
with the literature (see e.g. [131 US]) in a separate publication [48 . 

We can try to measure these numbers numerically using the method described in section [31 To 
do so, we need to construct the lattice version of the Aj operator using eq. (|3.6p with = e^, 
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L^2N 


/3l,5 


I3l.7 


P1..9 


8 


-0.937759 


-13.3574 




10 


-0.959708 


-14.8908 


-1518.37 


12 


-0.971844 


-15.7936 


-1805.43 


14 


-0.979236 


-16.3612 


-2013.66 


16 


-0.984064 


-16.7384 


-2157.86 


18 


-0.987388 


-17.0006 


-2262.59 


20 


-0.989771 


-17.1898 


-2340.51 


22 


-0.991539 


-17.3304 


-2399.80 


24 






-2445.81 


oo 


-1.0000 ± 0.0002 


-18.0(0) ± 0.05 


-27(00) ± 25 


Exact 


-1 


-18 


-2700 



Table 1: Measure of indecomposability parameters in the XXZ spin chain at q — i. 

hoo — 2/tt and vp — 2. For j — 2, 3, 4, these operators read 

A2 = i-i, /3i,5 = -l (4.8a) 
^3 = i^i-2L_2, /3i,7 = -18 (4.8b) 
Ai = L\ - 8L_2i-i + 12L-3, /9i,9 = -2700 (4.8c) 

We measured these three parameters in the spin sectors Sz = 1,2 and 3, respectively. The results are 
presented in Tab. [1] in very good agreement with the theoretical expectation. We computed these 
numbers using both supersymmetric (XX) and geometrical representations of the Temperley-Lieb 
algebra and obtained the same values in finite size. 

4.2 The cases x = 2 (c = 0) and x = 1/2 (c = —7): Percolation and st(2|l) 
superspin chain 

We now deal with a slightly more complicated example, the antiferromagnetic sl(2|l) superspin 
chain. This chain is known to be equivalent to the classical percolation problem, and arises naturally 
in the context of the Spin Quantum Hall Effect [1] . The percolation problem has c = and can 
also be formulated in terms of geometrical clusters or the XXZ spin chain at g = g2j7r/3^ -y^g ^^sq 
consider the ferromagnetic version, with central charge c — —7. 

4.2.1 s[(2|l)-invariant spin chain 

We consider a chain of alternating fundamental and dual representations of the Lie superalgebra 
s((2|l) with 'Hilbert' space "H = (□ (8) Q)®^. For more details about s[(2|l) and its representation 
theory, we refer the interested reader to the literature (see e.g. Ref. [13] )• The Hilbert space on 
one specific site is spanned by three independent states so the whole Hilbert space has dimension 
3L _ gjv^ every site, we introduce two boson operators [bi^c!,b'^j = SijSaa', where a S {t,i}, 

and one fcrmion {/i, /J} = (— l)*+^(5ij. We add a constraint to the system so that there cannot be 

more than one particle by site. The representation on the site i is C'^ ~ Span{// |0) , |0) , b\ ^ |0)} 
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Percolation {q ~ 


/3) 


Ferros[(2|l) [q^e^^-^l'i) 


L 




L 




L 




7 


-0.471874 


8 


-0.609088 


8 


-1.57616 


9 


-0.476386 


10 


-0.605858 


10 


-1.77278 


1 1 
ii 


-0.479983 


12 


c\ Cic\ci /I no 

-U.6U64UO 


12 


-1.87138 


13 


-0.482724 


14 


-0.607775 


14 


-1.92567 


15 


-0.484837 


16 


-0.609226 


16 


-1.95770 


17 


-0.486503 


18 


-0.610561 


18 


-1.97762 


19 


-0.487845 


20 


-0.611738 


20 


-1.99049 


21 


-0.488946 


22 


-0.612764 


22 


-1.99906 


oo 


-0.5000 ± 0.0001 


oo 


-0.6249 ± 0.0005 


oo 


-2.00 ± 0.005 


Exact 


-1/2 


Exact 


-5/8 = 0.625 


Exact 


-2 



Table 2: Measure of indecomposability parameters in c = and c — —1 theories. 

and corresponds to the fundamental □ for i odd and to the dual □ otherwise. The invariant couphng 
is chosen to be the Casimir in the tensor product representation of the sites i and « + 1 . 

= (^l+i,;&It + ^li.T^a + (-l)'+'/^i/^)(^.,Tfe^+i.i + Kxh^i.^ + (4.9) 
It furnishes a representation of the Temperley-Lieb algebra with n = I. The Hamiltonian reads 

L-l 

H = ±J2e^, (4.10) 

i=l 

where the minus sign corresponds to percolation (c = 0) and the plus sign to the ferromagnetic case 
(c = —7). There are two good quantum numbers Sz and B that we can use to label the states 

1 ^ 

= 2T.(^hf'^.^-^lM^^ (4.11a) 

1=1 

Finally, note that the low energy physics of this chain can be described by a non-linear sigma model 
with target space CP^'^ a,t 9 = tt [5T]. 

4.2.2 Measure of indecomposability parameters 

First of all, let us focus on the first known indecomposability parameter which concerns the stress 
energy tensor in the percolation problem. The Kac formula with x = 2 reads 

and the values appearing in the spectrum are 

/ii,i+2j = 5 • (4.13) 
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The partition function of the q = e"/^ XXZ spin chain with an even number of sites reads 

°° (4j-l)V24 _ (4j+5)V24 

Z = + I)- . (4.14) 

For the s[(2|l)-invariant chain, one replaces 2j + 1 by [2j + with q' + q'~^ = 3 [T51I32]. As we 

aheady argued in details in section 11.21 the stress energy tensor at c = must have a logarithmic 
partner t{z) which corresponds to the field ip{z) with our notations. It is now well admitted that 
the indccomposability parameter in this case is /3i.5 = —5/8. This number was also measured 
numerically in Ref. ^20j . The Jordan cells equations in that case read 

Lo \T) - 2 |T) (4.15) 
Lo \t) = 2\t) + \T) 

\T) = L-2 |0) 
L2\t)=b\0). 

Using the general structure (section (231) , we see that these states are organized into the following 
diamond structure 

i?2 t 
^ \ ^ \ 

V2= Ro Rs = I p (4.16) 

\ / \ ^ 

We measured /3i^5 for various Temperley-Lieb representations. In the XXZ spin chain ai q = e*'^/'^, 
the Jordan cell occurs in the sector Sz — 0, while in the sl(2|l) SUSY case, the Jordan cell for T is to 
be found in the sector (S'^, B) = (0, 0). We recall that the lattice indccomposability parameter given 
by eq. (|3.7p does not depend on the chosen representation. In the geometrical setup of percolation 
as dense loop gas with fugacity n — \, the Hamiltonian remains fully diagonalizable and there is 
no coefficient to measure here. Nevertheless, it is still possible to slightly deform it [5D] so that 
Jordan cells appear, and in this case we find the same values as in the other representations. The 
results (Tab. ^ are in excellent agreement with the prediction b = —5/8 and significantly improve 
the precision of the results obtained from the trousers trick [2D] . 

One can also measure the indecomposability parameters Pi ^1+2 j with j half-integer from odd- 
length chains. For instance, for L odd, there is a Jordan cell at level 1 that corresponds to the 
Virasoro staggered module V^. This Jordan cell occurs in the sector — 1/2, we call as usual ^ 
the unique state with /i = in this sector, and and </> the states with h — 1. We normalize ^ such 
that (flO = 1. The OPE formula gives a logarithmic couphng /3i_4 = {(t)\i}) = -1/2; this 

value was also found by Mathieu and Ridout using different methods [2i [13] . We can measure this 
coefficient on the lattice using the same method; once again, the results are in excellent agreement 
with the theoretical expectation (Tab. [2]). 

Another case of interest is g = e^"/"^, with a central charge c = —7. This model corresponds to 
an s[(2|l) spin chain with ferromagnetic couplings. Graphical expansion of the partition function 
yields a loop model with fugacity n = — 1. The conformal dimensions appearing in the spectrum 



-"^"We define the q-analog as [n], 
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L = 2N 


Pl,5 


8 


-1.26986 


10 


-1.29548 


12 


-1.31743 


14 


-1.33489 


16 


-1.34876 


18 


-1.35993 


20 


-1.36905 


22 


-1.37663 


oo 


-1.4582(8) ± 0.0001 


Exact 


-35/24 - -1.4583 



Table 3: Measure of (5 1^5 in the XXZ spin chain at q = e 



are given by the Kac formula for a; = 1 /2 



hi,i+2j — 



Jjj - 2) 
3 ' 



(4.17) 



The first interesting Jordan cell arises at level 1, and corresponds to the staggered module 



V3 



R2 



R5 



(4.18) 



Rs 



In this case, eq. (|1.22p with A = L_i yields = —2. This coefficient was also computed by 
Kausch and Gaberdiel pjj thanks to the Nahm-Gaberdiel-Kausch algorithm. This cell occurs in 
the sectors Sz — — 2,— ljj^2 of the XXZ spin chain. Measures of /3ij in all these sectors yield 
values in good agreement!^ with /3i_7 = —2 (see Tab. [5]). As in the other cases, the lattice values 
of /3i_7 do not depend on the chosen representation. 

4.3 The case a; = 3 (c = |): Logarithmic Ising model 

Let us consider the case q = e"/**, which corresponds to a central charge c — ^. It corresponds to 
a dense loop model with fugacity n = ^/2; this is of course the celebrated (logarithmic) Ising model 
CM{3, 4). The spectrum is given by the Kac formula at x = 3 

j(3j ~ 1) 



hi 



l+2j 



The partition function of the U„^^i^/4,{sl2) XXZ spin chain reads 



Z = ^(2j + 1 



(4.19) 



(4.20) 



-'^-'^ Extrapolation arc done fitting data by /3 + A/L + B/L^ , so the resulting fitting curves need not be monotonic; 
this is particularly obvious in this case. 
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L = 2N 




4 


0.021029 


6 


0.145101 


8 


0.276585 


10 


0.382046 


12 


0.463292 


14 


0.526436 


oo 


0.9 ± 0.1 


Exact 


5/6 ~ 0.8333 ... 



Table 4: Measure of fS^^i in the 0{n 0) model. The convergence is quite poor compared to the 
Temperley-Lieb case, but the result is not so bad given that the precision we obtain is roughly the 
same as for the first critical exponents. The result is consistent with the known value = §• 



The algebraic structure of the continuum limit was studied in details in section 12.31 We measured 
numerically the indecomposability parameter associated to the module (|2.29p . The results are shown 
in Tab.[H]and are in excellent agreement with the expected result /3i_5 — —35/24 (see eq. (|2.35p ). 

4.4 Another theory at c = 0: 0{n — j- 0) model and dilute polymers 

Finally, we present here an example of measure of indecomposability parameters in 'dilute' LCFTs. 
We study the 0{n — > 0) model, which is known to be relevant for the physics of dilute polymers. 
We can also formulate this model using supersymmetry, in terms of osp(2TO|2TO)-invariant spin 
chain and to non-linear sigma models with supersphere target space 5'2™-i|2m r^-j^Q partition 
function of the S ~ I C/g(s[2)-invariant chain reads 

o'»i+2j,l _ „'i-l-2j.l °° (l+6j)V24 _ (5+6j)V24 

where the exponents appearing in the spectrum are now in the first column hi^2j,i of the Kac table. 

We are interested in the stress energy tensor T in this theory, which is primary in this case 
as the central charge is c = 0. As in the percolation problem, T has a logarithmic partner that 
we call t. However, while in the percolation theory t had to be somehow identified with the field 
$i_5(z), it corresponds here to $3^1 (z) with thus a completely different indecomposability parameter 
Ps,! = 5/6. We would like to measure this coefficient from the lattice model defined eq. (|2.19p . 
Unfortunately, the finite-size convergence of the exponents is very poor, even the central charge 
cannot be properly measured in this case. This is probably due to the absence of integrable K- 
matrices in our model; a similar phenomenon was observed in Ref. [H]. Of course, one could add 
K-matrices at the edges in order to improve the convergence. Nevertheless, it is not clear to us 
how to adapt eq. (|3.6p in that case. Note also that the Hilbert space is much larger here so the 
accessible sizes are relatively small. 

Nevertheless, we can still hope to deduce a rough estimate of the indecomposability parameter 
b for the stress energy tensor here. Since the action of the operator L2 on the vacuum gives 0, we 
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define the operator yl^^^ as 



(L) 



^-2 — ^ 



2L 



t:vf 



{hi - hoo) cos 



2j7r 

~17 



(4.22) 



Of course, we could also have used A^^^ = but as it turns out, the convergence is better with 



-(L) 



+ L^^2- Measure of b from this formula arc shown in Fig. 21 Although the convergence is 



clearly not as good as for the previous examples, the result is consistent with the value 5/6. 



4.5 Remarks on descendants 

We have seen that using the structure of staggered modules over Virasoro which arise in a given 
theory, one can predict the whole structure of Jordan cells in the continuum limit. However, this 
does not take into account descendants: there is a whole pyramid of Jordan cells associated with 
each (primary) Jordan cell. To be more precise, if there is a Jordan cell mixing two operators 
if) and (j) with parameter /3, then we can expect Jordan cells for all the descendants as is readily 
shown using the commutation relations of the LnS. The resulting indecomposability parameters 
are not independent and can be deduced from the knowledge of /3. For instance, let us consider the 
operators V'*-""'' = ^-nV' and = L_„0 in the case A = L_2- In the basis (■!/''■""', 'Z'^""''), Lq 

reads 

where /i is the conformal weight of -0 and (/). Let /3^~") = <^?/;(~")|(/)(^")^. Using the Virasoro algebra, 
one can show that this coupling is given by 

= {-^^n^ - 1) + l^na + 4n) p. (4.24) 

Of course, there are similar formulae for other kinds of descendants and for other A operators. 
We remark that these results are compatible with eq. (|1.22p . We measured indecomposability 
parameters for descendants in some cases (results not shown here) and found a good agreement 
with the previous considerations. 



5 Conclusion 

Pushing further the analysis of [T^, we have shown in this paper that it is possible to investigate 
the fine structure of indecomposable Virasoro modules in LCFTs using numerical analysis of a 
certain type of lattice models. Our method is general enough to be adapted to many cases, and the 
precision reached is almost as good as for critical exponents. 

We have restricted to the simplest type of boundary conditions for the LCFTs, but extension to 
more complicated cases is possible using more complicated lattice models, based for instance on the 
blob algebra [13113117] (or the one and two-boundary Temperley-Lieb algebra). More interestingly 
maybe, we believe that the method can be extended to the periodic case as well and thus should 
provide a powerful tool to investigate the structure of bulk LCFTs, where very little seems to be 
known at present. We will report on all these questions soon. 
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35 
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693 
100 


6114399291 
1078465600 
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3-state Potts £M(5,6) 
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Dilute LCFTs 
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5 
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67375 
676 
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106462606250 
116550867 
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141745038705442046875 
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Going back to the values of the indecomposabihty parameters, we also argued that they can be 
inferred from a simple heuristic argument relying on OPEs. This did not seem to be known, and 
suggests revisiting the bulk problem as well, by systematically considering LCFTs as the limit of 
usual, non-logarithmic, CFTs. This will also be discussed elsewhere. 

Finally, we summarize our results with a table of the first few indecomposabihty parameters for 
the minimal logarithmic models CM{p,p'). We focus on the series CM{l,p) and £M{p,p + 1), 
and we consider the two versions 'dense' and 'dilute' of each theory. Using the general structure 
of section 12.41 and eq. (|1.22|) , one can analyze the indecomposabihty parameters in a systematic 
fashion [IS]. The operators A = L_n + ■ ■ ■ are generated using the null- vector condition, and are 
normalized as in eq. (|1.4p . The results are gathered in Tab. 

For a given dense {resp. dilute) CM{p,p') theory, we denote (3i,i+2j {resp. /3i+2j,i) the log- 
arithmic coupling associated with the Jordan cell at level /ii_i+2j {resp. /ii+2j,i). Of course, it 
may happen that this Jordan cell does not exist, or there may not be any interesting coefficient to 
measure (this is the case for the first few Jordan cells in £A4{l,p) theories), in which cases we use 
different symbols. Some of the parameters we find are quite complicated irreducible fractions, and 
the simplicity of the results sometimes depends on the normalization choice for A. Note that this 
table contains only a small fraction of the couplings that eq. (|1.22p allows to compute. In principle, 
the formula p.22p could be applied to obtain any indecomposabihty parameter for a given theory. 
The limitation obviously comes from computing Virasoro commutators. We shall report on all this 
in gg. 
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